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Q ■ ABSTRACT 

\ We present Herschel-PACS spectroscopy of the [OI]63 /im far-infrared cooling line from 

a sample of six unlensed and spectroscopically-confirmed 870 /im-selected submillimetre 
(submm) galaxies (SMGs) at 1.1 < z < 1.6 from the LABOCA Extended Chandra Deep 
Field South (ECDFS) Submm Survey (LESS). This is the first survey of [OI]63 /im, one 
■ of the main photodissociation region (PDR) cooling lines, in SMGs. New high-resolution 

\^ ' ALMA interferometric 870 fim continuum imaging confirms that these six Herschel-taigeted 

, SMG counterparts are bona fide sources of submm emission. We detect [OI]63 /im in two 

SMGs with a SNR> 3, tentatively detect [OI]63 fim in one SMG, and constrain the fine flux 
for the non-detections. We also exploit the combination of submm continuum photometry 
from 250-870 fLm and our new PACS continuum measurements to constrain the far-infrared 
luminosity, Lfir, in these SMGs to < 30 per cent. We find that SMGs do not show a deficit in 
their [OI]63 /im-to-far-infrared (FIR) continuum luminosity ratios (with ratios ranging from 
0.5-1.5 per cent), similar to what was seen previously for the [CII]158 /im-to-FIR ratios 
in SMGs. These observed ratios are about an order of magnitude higher than what is seen typ- 
ically for local ultra luminous infrared galaxies (ULIRGs), which adds to the growing body 
of evidence that SMGs are not simply 'scaled up' versions of local ULIRGs. Rather, the PDR 
line-to-LpiR ratios suggest that the star formation modes of SMGs are likely more akin to 
that of local normal (lower-luminosity) star-forming galaxies, with the bulk of the star forma- 
tion occurring in extended galaxy-scale (^kpc) regions. These observations represent the first 
step towards a census of the major PDR cooling lines in typical SMGs that will be attainable 
with ALMA, enabling detailed modelling to probe the global properties of the star formation 
and the evolutionary status of SMGs. 

Key words: galaxies: high-redshift - submillimetre: galaxies - galaxies: star formation - 
galaxies: starburst - galaxies: ISM - infrared: ISM 
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*Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic- 
ipation from NASA. 



1 INTRODUCTION 

Far-infrared bright galaxies comprise some of the most luminous 
galaxies in the local Universe, but they contrib ute just 0.3 per cent 
of th e total local galaxy luminosity density jSanders & Mirabell 
Il996h . At first glance it appears that they are an unimportant com- 
ponent of the galactic zoo, however, sensitive (sub)millimetre cam- 
eras have revealed that ultraluminous infrared galaxies (ULIRGs; 
Lfir > 10^2 L©) are much mor e abundant at high redshifts 
{z > 2). Since their discovery {Smail, Ivison & BlainI 1 19971 : 
iHughesetalll 19981 ; ISarger et alj|l998l) . nearly 13 years of piece- 
meal survey work has so far produced samples tot alling ^ 500 
850-12Q 0-selected (sub)mm galaxies (S MGs; e.g. ICoppin et al.l 
[2006: Austermann et alJlOlOl l lWeifi et al.ll2009.) . wit h extreme star 
forma tion rates (SFR s) of ~ 100-1000 yr~^ ( e.g. lXovacs et al.l 
l2006l : ICoppin et al.l l2008l : iMagnelli et all I2OI2I) and fuelled by 
vast reservoirs of cold molecular ga s (Mh, 5 X 10^° Mr ; 
iGreve et alj2005l:lTacconi et al .'2006. 2008: B othwell et al.ll2012t) . 
Spectroscopic follow-up of a subset of SMGs suggests that the vol- 
ume density of ULIRGs increases ~ 100 times out to z ~ 2 
jChapman et alj l2005l : iBanerii et al.l 1201 ll) . Therefore, luminous 
obscured galaxies likely dominate the total bolometric emission 



from star formation at early epochs (e.g. IPascale et ai] I2OO9I : 
lMarsdenetai]|2009l : lEales et alj|2009l . I2OIOI : lOliver et al.l l2010l) . 
and possibly correspond to the a ctive starburst pro genitor phase of 
massive spheroidal galaxies (e.g. lLillv et al Il99 9^. 

However, several critical questions about SMGs remain unan- 
swered. The most fundamental of these is: How is the immense 
luminosity of an SMG powered? - by an intense single compact 
nuclear starburst (as seen in local ULIRGs) or by a more ex- 
tended starburst (as seen in local normal star-forming galaxies)? 
One route to tackle this question is to ask what are the local ana- 
logues of SMGs and can these be used to understand the physics 
of the star formation activity within the distant SMG population? 
Are SMGs simply 'scaled up' ULIRGs, where the star formation 
occurs in a compact and highly obscured nuclear region - poten- 
tially with a significant contribution from an active galactic nucleus 
(AGN; e.g. [Alexander et alj|2005bl) ? Or does some fraction of the 
population resemble very scaled up 'normal' star-forming infrared 
(IR) galaxies, which have far-IR em ission arising from an extended, 
cooler component (termed 'cirrus' : [Efstathiou & Rowan-Robinsonl 
I2OO3I) ? If we could empirically match subsets of the SMG popula- 
tion to local galaxies, then this would provide significant insights 
into the physical processes occuring within this important popula- 
tion. The synergy of Herschel and the Atacama Large Millimetre 
Array (ALMA) provides a novel opportunity to move beyond the 
first generation of detection experiments by beginning to examine 
the interplay between the interstellar medium (ISM) and the power 
source of high-redshift luminous dusty star-forming galaxies. 



1.1 Far-infrared spectroscopy: a probe of SMG astrophysics 

The primary diagnostic lines of the ISM in the far-IR have tra- 
ditionally included the fine-structure lines of [OIII]52, [NIII]57, 
[OI]63, [OIII]88, [NII]122, [OI]146, [CII]158^m, [CI]370Atm, 



and [CI]610 ^m. The relative strengths of these lines can be mod- 
eled as a function of the neutral gas cloud density and the cloud- 
illuminating far-UV radiation field intensity in p h otodissociation 
regions (PDRs; e.g. Tielens & Hollenbach 1985: Kaufm an et"al] 
ll999l : IWolfire et al.l2003 : Meijerink & Spa ans 2005), providing di- 
agnostics of the physical conditions of the ISM. While much has 
been learned about the ISM conditions in Galactic star-forming 
regions and local galaxies th rough facilities such as the Infrared 
Space Obse r vatory (ISO; e.g. iMalhotra et al]|200ll : iLuhman et al.l 
1 19981 I2OO3I : iNegishi et all I2OOII) and the K i iiper Airborne Ob- 
servatory (KAO; e.g . ICrawford et al.l 1 19851 : Istacev et alj Il99ll: 
IPoglitsch et aljri99a) , due to previously poor instrument sensitiv- 
ity and inaccessibility from the ground, dedicated surveys prob- 
ing the ISM conditions in samples of high re dshift dusty galax- 
ies have only rece ntly begun in earnest (e.g. IStacev etal.ll201(]|: 
iDecarU et ani2012h with the commissioning of purpose-built in- 
strumentation like the redshift (z) and Early Universe Spectrome- 
ter (ZEUS; Stacev et al. 2007; Hailev-Dunsheath 2009). For galac- 
tic and extragalactic systems, the 'classical' PDR tracers are the 
brightest fine-structure emission lines [CII]158 pm and [OI]63 pm, 
acting as the primary coolants for the dense (n > 10-10^ cm~'^), 
warm (T~ 100-1000 K), and neutral material. Since [CII]158^m 
and [OI]63/im correlate well with SFR and are sufficiently lumi- 
nous (accounting for as much as ~ 0.1-1% of the bolometric lu- 
minosities in nearby gala xies and are th us much more luminous 
than molecular lines; Sta cev et"al]|l99lh . they have the potential 
to be used as powerful tracers of ISM conditions out to high- 
z. However, ISO observations of local ULIRGs revealed that the 
[CII]158 ^m emission is about an order of magnitude lower rela- 
tive to the far-IR (FIR) continuum flux than is seen in normal an d 
starburst galaxies (e.g. IMalhotra et alj|2001 ; Luhman et al.ll2003h . 
And this has recently been shown by Graciii-Carpio et alj ( 1201 ih 
using Herschel-FACS to be a general aspect of far-IR fine struc- 
ture lines (see Section l4l and iGracia-Carpio et alj|201 ll for a dis- 
cussion of a possible physical explanation). The potential weak- 
ness of [CII] 158 fim has implications for its use as a star formation 
tracer for hi gh-redshift studies (z > 4; e.g. Maiolino et al.ll20o3 : 
although cf. lMaiolino et alj|2009l : iDe Breuck et alj|201 ih . In order 
to make further progress in understanding the driver of star forma- 
tion in luminous high-redshift SMGs in relation to galaxy popu- 
lations in the local Universe, we need to begin compiling a simi- 
lar set of (multi-transition) CO, [CII] 158 ^m, and [OI]63 /im data 
(which provide strong joint constraints on the PDR and star for- 
mation conditions such as the Hydrogen volume density, n, and 
the far-UV radiation field strength. Go in Habing units) for typical 
luminous star-forming SMGs at the epoch where their population 
peaks (z ^--^1-3). 



1.2 This paper: New Herschel observations of [OI] 63 /im in 
SMGs 

In the present paper we will be discussing observations perform ed 
with the ESA Herschel Space Observatory jPilbratt et ani2010l) . in 
particular employing Herschel 's large telescope and powerful sci- 
ence payload to perform photometry and spectrosco py using the 
Spect ral and Photometric Imaging REceiver (SPIRE: [Griffin et aU 
l2010t) and the Photodetecto r Array Camera and Spectrometer 
(PACS:'Poglitsc h et al.ll2010 Y This paper presents the first attempt 
at compiling a census of the strength of [01] 63 jim in far-IR lu- 
minous SMGs at the peak of star formation activity in the Universe 
(z ~ 1-2), enabling us to determine its suitability as a star formation 
tracer at high redshift. These Herschel observations represent a key 
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stepping stone to future ALMA studies of [CII] 158 /im and (multi- 
transition) CO to investigate the ISM physics of SMGs in detail, 
since [OI]63 /im emission can only be measured in 1 < z < 2 
SMGs with Herschel-FACS, as the emission line falls outside of 
the ALMA bands for all but the highest redshifts {z > 4). 

This paper is organized as follows: the sample selection, 
Herschel-FACS spectroscopic observations and data reduction are 
described in Section[2] In Section[3] we present the main analysis 
and results of the PACS spectrocopy and combine these measure- 
ments with Lfir estimates obtained f rom the combination of Large 
Apex Bolometer Camera (LABOCAi lSiringo et al.ll200si) and new 
Herschel-SPIRE and PACS photometry. We discuss the implica- 
tions of our results in Section|4] Finally, our conclusions are given 
in Section [5] Throughout the paper we assume cosmological pa- 
rameters of Q,a = .73, = 0.27, and _ffo = 71kms"^Mpc"^ 
dSpergel et alj2003l) . 



ECDFS/CDFS iLehmer et al.l 120051 : IXue et alj 1201 ih . Here we 
also make use of new public Herschel Multi-tiered E xtragalactic 
Survey (HerMES; iHerMES Collaboration et al. Illoil Herschel- 
SPIRE maps to extract 250-500 /im photometry in combination 
with the zLESS spectroscopic redshifts, providing the tightest con- 
straints on Lfir currently possible for these SMGs (to < 30 per 
cent). In addition, high-resolution (~ 1.55" full width at half maxi- 
mum; FWHM) ALMA 870 ^m continuum data have been collected 
for the majority of the 126 LESS SMGs during cycle to a 1 a rms 
depth of ~0.4 miy, yielding positions accurate to < 0.3" (Karim et 
al., in preparation). The ALMA data have enabled us to verify un- 
ambiguously the validity of the radio and/or mid-IR counterparts 
for our SMG sample, which are used here as a basis for the PACS 
target positions and central frequency tuning of the spectrometer 
(see Section |3.1.1| for a summary of the ALMA ID results for our 
sub-sample). 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Sample selection and multi-wavelength properties 

The requirement of our survey is a sufficient number of SMGs se- 
lected in a uniform manner that have been spectroscopically con- 
firmed to lie between 0.7 < z < 2 so that the redshifted [OI]63^m 
emission line will lie within the most sensitiv e region of the first or - 
der of the PACS spectrometer (107-189 Mm JPoglitsch et alj2O10h . 
The SMGs also need to be accessible by ALMA for future stud- 
ies (i.e. S < 30°). One such field is that covered by the 870 ^m 
LABOCA Extended Chandra Deep Field South (ECDFS) Submil- 
limetre Survey (LESS; WeiB et al. 2009), which comprises a parent 
sample of 126 SMGs detected above 3J-a in the largest contigu- 
ous (30' X 30') a nd uniform map with a 1 a rms of ~ 1.2 mjy. 
lBiggsetal.l ( l201ll) have identified secure radio and mid-IR coun- 
terparts for 63 per cent of these SMGs. Since the ECDFS is one 
of the pre-eminent cosmological deep fields and has as such been 
extensively imaged and targeted spectroscopically, some of these 
counterparts already possess secure archival redshifts. Redshifts for 
the remaining counterparts without archival spectroscopy were tar- 
geted through a large observational programme (zLESS; Danielson 
et al., in preparation) using the European Southern Observatories 
(ESO) Very Large Telescope (VLT) FOcal Reducer and low dis- 
persion Spectrograph v.2 (F0RS2) and Visible MultiObject Spec- 
trograph (VIMOS). At the time the Herschel proposal was submit- 
ted in 2010, the VLT survey was ~ 40 per cent complete. The 
total number of SMGs fulfilling the requirements outlined above 
at that time was nine, and these nine SMGs thus comprise our 
Herschel-PACS spectroscopy sample. Thus, our sample is com- 
plete in terms that it is comprised of all of the spectroscopically 
confirmed 1 < 2: < 2 SMGs in LESS that are radio and/or 24 ^m 
identified. 

We make use of additional multi-wavelength data sets avail- 
able for the radio and/or 24 fim LESS SMG counterparts in or- 
der to further interpre t our results, which we now briefly describe. 
IWardlow et alj ilOl ih measured the optical-mid-IR photometry for 
the robustly identified LESS SMG counterparts and provide photo- 
metric redshift estimates and Lir estimates. In order to assess any 
potential AGN contribution i n the SMGs , we use the Spitzer In- 
fraRed Array Camera (IRAC; Fazio et al. 2004) photometry from 
IWardlow et all ilOUh extracted from the Spit zer IRAC/MUSYC 
Public Legacy Survey in the ECDFS (SIMPLE; lDamen et al.l201lh 
and the 250 ks and 4 Ms X-ray imaging and catalogues of the 



2.2 Herschel-PACS spectroscopy observations 

The [0I]63 /im emission line in our sample of SMGs was observed 
with the Herschel Space Observatory' PACS red camera over 112- 
167 /xm using the high-resolution range scan spectroscopy mode 
and employing a chop/nod to provide background subtraction. For 
each target, we covered a small spectral range of ~ 4/im, not- 
ing that only the central ~ 3 /xm of the spectra are usable, corre- 
sponding to A2 ~ ±0.03-0.05 to account for small redshift errors 
and/or systematic offsets from the optical redshift, as well as to al- 
low for a reliable continuum or baseline fit. The aim was to reach 
a uniform 1 a line sensitivity of ~ 0.8 x 10~^* Wm~^ across the 
sample (see Table reftabiobs for the actual sensitivities that were 
reached). Six targets were observed in 2012 February over a total of 
15.4hrs (including overheads). Table[T]gives the target list, obser- 
vation mode parameters, and individual integration times. We note 
that an additional three sources were targeted as part of this pro - 
gramme (LESSIO, LESS50 and LESS118; see lBiggs et a"i]|2010l) . 
but the IDs (and therefore redshifts) have been subsequently revised 
based on new ALMA continuum data (Karim et al., in preparation; 
see also Section [3.1.1| l. Thus, these three observations do not pro- 
vide interesting constraints on the [0I]63 /xm emission from SMGs 
and so we have dropped them from the analysis. 

2.3 Data reduction 

The data were reduced w ith the H erschel interactive processing en- 
vironment (HIPE v8.1.0; lot^l201Q) . We employed the background 
normalisation pipeline script distributed with HIPE, which imple- 
ments an alternative way to correct for drifts in the response of 
PACS during the observation and an alternate flux calibration than 
the standard pipeline. Using this alternative pipeline script results 
in marginally better quality data products with smoother continua, 
which is helpful when searching for faint spectral features atop a 
faint continuum. As part of the standard pipeline, the continuum 
in each of the 16 spectral pixels was scaled to the median value 
to correct for residual flat field effects, and finally the two nod po- 
sitions were combined to completely remove the sky (telescope) 
background. Given that our targets are point sources to PACS we 
have measured line fluxes from spectra extracted from the central 
9.4" X 9.4" spatial pixel of the 5 x 5 pixel field-of-view of the PACS 
spectrometer. The pipeline automatically applies a beam size cor- 
rection factor and an additional in-flight correction of the absolute 
response. The spectra have been re-binned to Nyquist sample the 
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Table 1. Observation details and measured [OI]63 integrated lin e fluxe s and continuum flux densities from our Herschel-¥ACS programme 
(OTl_kcoppin_l). SMG optical counterpart positions from lWardlow et al] (2oTl|) are given in columns 3^, to which the zLESS and archival spectroscopy 
are referenced. The op tical spectroscopic redshifts in column 5 are from: (a) zLESS (Danielson et al., in preparation); (b) Vanzella et al. (2005); and (c) 
ISilverman e t al. ( 2010). Columns 6-9 give the Herschel Observation Date (OD), Observation identification number (ObsID), the number of range repetitions 
per nod cycle times the number of nod cycles, and the total Astronomical Observation Request (AOR) duration (including the slew). Column 10 gives the mea- 
sured (continuum-subtracted) integrated [OI]63 fj,m flux and error from the Gaussian fit described in Section lsTT] or the 3-cr upper limit assuming a Gaussian 
profile with a FWHM= 300 km s ^ ^ and a peak height of 3 X crrms . Column 1 1 gives the best-fitting continuum flux density (or upper limit) at the spectrum 
centre (indicated in square brackets). Column 12 gives the 1 a rms of the entire spectrum, quoted for the representative bin width at the centre of each band 
(see square brackets in previous column). We note that the line and continuum flux density error budgets do not include the estimated additional absolute 
calibration rms uncertainty of 12 per cent (see text). 
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LESSJ033329.9- 


-273441 


LESS21 


3:33:29.73 


-27:34:44.4 


1.235° 


1020 


1342239701 


4 


X 


2 


8146 


< 3.8 


< 18 [141 /jm] 


3.9 


LESSJ0332I7.6- 


-275230 


LESS34 


3:32:17.61 


-27:52:28.1 


1.0979*' 


1020 


1342239703 


4 


X 


2 


7874 


< 3.9 


< 17 [133 urn] 


4.1 


LESSJ033331.7- 


-275406 


LESS66 


3:33:31.92 


-27:54:10.3 


1.315^^ 


1013 


1342239369 


4 


X 


2 


8349 


2.0 ± 0.6 


27 ± 7 [146 Mm] 


4.1 


LESSJ033155.2- 


-275345 


LESS88 


3:31:54.81 


-27:53:41.0 


1.269" 


1020 


1342239705 


4 


X 


2 


8214 


1.8 ± 0.6 


23 ± 7 [143 Mm] 


4.2 


LESSJ033140.I- 


-275631 


LESS106 


3:31:40.17 


-27:56:224 


1.617" 


1019 


1342239753 


5 


X 


2 


11444 


3.8 ± 1.4 


< 23 [165 fim] 


44 


LESSJ033150.8- 


-274438 


LESS 114 


3:31:51.09 


-27:44:37.0 


1.606" 


1020 


1342239702 


5 


X 


2 


11444 


< 5.6 


48 ± 11 [165 pm] 


5.3 



spectral resolution element, which is appropriate for resolved lines 
(see Fig.[Tll. We have chosen not to upsampl^the data, so that each 
bin is independent, yielding reliable estimates of the noise proper- 
ties in the spectra. A similar data reduction strategy was adopted by 
the Survey with Herschel of the ISM in Nearby Infrared Galaxies 
(SHINING) team for reducing faint spectra of high-redshift galax- 
ies (priv. comm. E. Sturm and J. Gracia-Carpio). 

Third order data in the wavelength range 51-73 fim are col- 
lected by the blue camera in parallel with the first order [01] 63 fim 
data. However, these observations are not frequency-tunable and 
cover only ~ 2 /im. Nevertheless, the data were reduced using the 
same procedure as above and examined for any serendipitous line 
detections, but the useful parts of the spectra do not cover any use- 
ful known lines at our source redshifts and so we do not discuss 
these data further. 



3 ANALYSIS AND RESULTS 

We present new [OI]63 ^m data for the SMGs, including two new 
detections of the emission line, and one marginal detection. For 
the non-detections, upper limits are only useful if both the counter- 
part and redshift are correct: this Herschel spectroscopy program 
was designed to target only the most statistically probable 'robust' 
counterparts also posessing good quality redshifts. However, given 
the significant time lag from the survey definition to the time of data 
acquisition (~ 1.5 years), new data have since been accumulated 
from zLESS and ALMA that can help us to verify if the targeted 
counterparts and redshifts originally chosen are still correct. In the 
following section we examine and discuss each PACS spectrum. 
We then derive the FIR continuum flux and Lfir for each SMG via 
a modified blackbody spectral energy distribution (SED) fit to pub- 
licly available submm data so that we can explore the [OI]63 fj,m- 
to-FIR continuum ratios of our sample as a function of Lfir in 
context with other low- and high-redshift galaxy populations. 



^ Upsampling refers to increasing the sampling rate of a signal. Here we 
use upsample=l. If upsample> 1 is used, then the frequency bins are not 
independent (ie. they will have some data in common) and this will have the 
effect of changing the spectral noise properties as well as the line widths. 



3.1 Spectral line and continuum measurements 

We detect weak [OI]63 fim emission in LESS66, LESS 106 and ten- 
tatively in LESS88 (see Fig.[B. We can derive useful upper limits 
for the cases where the line emission is not detected, and only in 
those cases where we are certain that the correct redshift and/or 
SMG counterpart was targeted. We detect continuum emission in 
LESS66, LESS88, and LESS 1 14 at > 3 cr, and for the continuum 
non-detections we derive 3cr upper limits to the continuum flux 
density in the PACS spectral range. 

For the detections, a Gaussian profile and continuum are fit- 
ted simultaneously to the Nyquist sampled spectra by performing 
non-linear Levenberg-Mar quardt least-square s minimization using 
the MPFIT IDL package l lMarkwardlll2009 '). Since the pipeline- 
calculated errors in the PACS spectrum are known to be unre- 
liable, we have clipped off the noisy edges from each spectrum 
(« 1 fim on each side) and weight each of the remaining central 
channel data points equally. To calculate the velocity-integrated 
flux of the [OI]63 /im emission line, we integrate the best-fitting 
(continuum-subtracted) Gaussian profile. The I a errors on the 
velocity-integrated flux and continuum flux density are derived 
from a bootstrap estimate generated by adding random noise with a 
I a rms equivalent to the I a rms of th e line-free channels, re peated 
1000 times with replacement (e.g. see Wall & Jenkins 2003). This 
should provide a reasonable estimate of the uncertainties since our 
choice of spectral binning (upsample=l) ensures that the noise 
should not be strongly correlated between channels. 

For the non-detections we calculate 3 a upper limits assum- 
ing a Gaussian profile with a 3 arms peak height centred at the 
expected position of the [OI]63/im emission line. We use the 
channel-to-channel I a rms noise across the clipped spectra and 
adopt an instrumental resolution-corrected emission line FWHM 
of 300 km s^^, which is consistent with our detections and those 
of ot her high-redshift g alaxies that have been observed with PACS 
(e.g. iSturmetalj20I(i) . 

Integrated emission line fluxes, upper limits and continuum 
flux densities are summarized in Table [T] Note that the rms uncer- 
tainties associated with the absolute calibration of PACS have not 
been folded into the quoted flux uncertainties and are estimated to 
be ~ 12 per cen{f|. 

^ See |http://herschel.esac.esa.int/twiki /bin/view/Public/PacsCaUbrationWeb| 
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Figure 1. Observed infrared Herschel-PACS Nyquist sampled spectra of our sample of six SMGs. The spectra have been zoomed to a 1 /im window centred 
on the optical or [OI]63 fim derived redshift, and the dot-dashed vertical line and eiTor bar indicate the expected position of the [OI]63 fj-m emission line 
based on prior optical spectroscopy and an approximate 1 cr error of 100 kms^^. For LESS66, LESS106 (formal detections) and LESS88 (only marginally 
detected) the best-fitting Gaussian+continuum profile is overplotted as a dashed curve (see Section lsHl for details). For the non-detections (LESS21, LESS34, 
and LESS 114) the best underlying continuum fit is overplotted as a dashed line. 



In the following we first summarize the counterpart identifica- 
tions in light of the new ALMA data, and then discuss the details of 
the counterpart identifications and redshift robustness for each tar- 
geted SMG in turn, along with the new PACS spectroscopy results. 
Reference to a 'robust' counterpart is meant in a statistical sense: 
given a potential radio or mid-IR counterpart with its flux density 
and radius from the SMG position, the a priori probability, p, of 
finding at least one object within that radius of at least that flux den- 
sity from the expected number of events at random is p < 0.05. A 



'tentative' counterpart refers to sources with 0.05 < p < 0.1 and is 
often still low enough to indicate a likely counterpart. Ilvison et"ail 
( 12002, , . 200?) and Biggs et al. ( 201 1) provide a more in depth expla- 
nation on this subject. Even though this is the only feasible route 
for precisely pin-pointing the positions of the bulk of SMGs in the 
absence of very high-resolution interferometric continuum submm 
data (ie. in the pre-ALMA era), this is only a statistical identifica- 
tion method. Recent submm inteferometric observations of SMGs 
have demonstrated that potentially a large fraction of 'robust' SMG 



© 0000 RAS, MNRAS 000, 000^000 



6 Coppin et al. 



identifications may be more complicated (with up to ^ 30 per cent 
breaking up into multiple components; se e e.g. ICowie et al] |2009l : 
IWang et al.ll20Tll : iHatsukade et"d]|201(]|) . Armed with new deep 
high-resolution ALMA 870 /im continuum data (Karim et al., in 
preparation), we can test the reliability of the counterpart identifica- 
tions using the radio or mid-IR, which is used here as a basis for tar- 
geting the [OI]63 yLva emission line with Herschel-FACS. Redshifts 
for the counterparts were obtained by an ESO VLT F0RS2A^IM0S 
Large Programme (zLESS; Danielson et al., in preparation) and 
are supplemented by the wealth of existing public archival spec- 
troscopy in the ECDFS. The quality of the redshifts are categorized 
by zLESS as follows: Q=l refers to a secure redshift (e.g. based 
on the identification of multiple lines) and Q=2 refers to a probable 
redshift (e.g. based on a single bright emission line identification). 

3.1.1 Summary of the counterpart ID re-examination in light of 
new ALMA data 

In summary, after examining the preliminary new ALMA contin- 
uum data from Karim et al. (in preparation), it appears that our 
original sample of nine targets (see Section l2Tt should be reduced 
to six. The ALMA data indicate that we have not targeted a coun- 
terpart associated with submm continuum emission in two cases 
(LESSIO and LESS118) and we thus discard these sources from 
our sample and from further discussion. The ALMA data have also 
not been able to validate the ID for LESS50 (whether this is because 
the submm emission is resolved out in the ALMA map or that the 
SMG has split up into several fainter sources when seen by ALMA 
- none of which correspond to the ID), and so to be conservative 
we also discard this source from the sample. Although, we note that 
our discarded targets are still radio or 24 /im emitting galaxies (just 
not the source of the submm emission) with good quality redshifts 
and thus remain useful data for other purposes. The six remaining 
targets are associated with ALMA continuum emission and pos- 
sess a good quality redshift so we can safely measure line fluxes or 
sensitive upper limits in these cases: LESS21, LESS34, LESS66, 
LESS88, LESS 106, and LESS 114. 

In a few cases (LESS34, LESS106, LESS114) it appears that 
the LABOCA-detected SMG may be splitting up into two or more 
continuum sources when seen by ALMA at much higher resolution. 
For these, ideally the LABOCA and SPIRE photometry should first 
be deblended or deconvolved before calculating the Lfir values in 
Section IT2I Unfortunately, without complete redshift information 
and multiple interferometric photometry for all of the ALMA con- 
tinuum sources this is impossible to do in a robust way, and it is 
beyond the scope of this paper. The diagnostic we are interested 
in comparing with literature values is the [01] 63 /im-to-FIR con- 
tinuum ratio versus Lfir for the SMGs. If the FIR continuum is 
reduced accordingly, the [01] 63 /im-to-FIR continuum ratio will 
increase, and so this has little bearing on our main finding in this 
paper that this ratio is high for SMGs compared to local galaxies 
with similar Lfir (see Section [373] l. Thus, by not correcting the FIR 
continuum luminosity we follow a conservative approach. 

3.1.2 LESS21 

It has a tentative 24 /im counterpart (e.g. borderline secure, with 
p = 0.067), and there is ALMA continuum emission coincident 
with this position, confirming the ID. 

The zLESS spectroscopic redshift of this counterpart is 2; = 
1.2350. The zLESS VIMOS spectrum of this counterpart is consis- 
tent with the photometric redshift solution oi z = 1.26lQ'5g3 frorn 



IWardlow et alj ilOl ll) . with the tentative detection of Mgll absorp- 
tion. It is not detected in the X-ray, and has a 3 o" upper limit on its 
full-band flux of So. 5-8 kcv < 6.1 x 10~^^ e rgcm~^ s~^, derive d 
from the 250 ks X-ray imaging of the ECDFS ("L ehmer et al.l2005l) . 
The S8;jm/S4.5^m colour ratio is ~ 0.7, indicating that the galaxy 
cont ains a negligible c ontribution to its energetics from an AGN 
(see lCoppinetalfcOlOl) . 

Given the optical redshift, the [OI]63 pm emission line is ex- 
pected to lie at 141.21 /im, however no line emission near this 
wavelength is apparent in the spectrum. We thus calculate a 3 cr 
upper limit to the emission line flux. 

3.1.3 LESS 3 4 

It neighbours LESSIO (~ 20arcsec away) and has a single se- 
cure 24 counterpart (and formally no secure radio counter- 
part), which was targeted with PACS. There does appear to be 
some weak ALMA continuum emission coincident with this tar- 
geted ID, however there is also another bright ALMA source lying 
just ou tside of the SMG counterpart radius adopted bv lBiggs et al.l 
( I20T l[). The arc hival redshift of our targeted source, z = 1.0979 
(Q=l), is from Vanzella et alj (|200^. The single SMG identified 
with LABOCA is thus possibly comprised of multiple sources of 
submm emission, which would have the effect of a reduced in- 
ferred Lfir for our targeted counterpart. It is not detected in the 
X-ray, and has a 3 upper limit to its full-band flux of So. 5-8 kcV < 
7.0 X 10~^^er gcm~^ s~^, de rived from the 4Ms X-ray imaging 
of the CDFS tXue et al.ll2oTTI) . The S8Mm/S4.5^,m colour ratio is 
~ 0.2, indicating that the galaxy c ontains a negligible contribution 
to its energetics from an AGN (see^Coppi n et alj201(t) . 

Given the redshift, we expect the [OI]63 nm emission line to 
lie at 132.55 /im. We do not see an emission line in the spectrum 
near this wavelength and so a 3 a upper limit to the emission line 
flux is calculated. 

3.1.4 LESS66 

It has a single robust radio and 24 /im counterpart. The ALMA con- 
tinuum data reveal continuum emission coincident with this ID. 
The zLESS redshift of this counterpart from an excellent qual- 
ity (Q=l) VIMOS spectrum is indicated to be 2 = 1.310, al- 
though zLESS F0RS2 spectroscopy indicates z — 1.314 (Q=l). 
A secure (Q=l) Keck-based spectroscopic redshift of z — 1.315 
was obtained by Silverman et al. (2010) from a program to tar- 
get optical and near-I R counterparts of faint X-ray sources in the 
ECDFS. Interestingly. ISilverman et ah I I2OIO.) classify this source 
as a Broad-Line AGN (BLAGN) based on the presence of at least 
one emission line having a FWHM> 2000 km s~^. It is detected 
in the X-ray with a full-band flux of So.s -skeV = (31.6 ± 1.4) x 
10~ ''' erg cm' s"^ iLehmeretaLllioC^ . Its X-ray flux is consis- 
tent with that expected for an optically classified QSO given its 
submm flux density and is hig her than typically seen in SMGs (~ 
lO"^*^ - 10-^"^ ergcm~^ s~ MAlexander et alj[20059 ). Although 
there is clearly a luminous AGN in this SMG, its S8^m/S4.5^m 
colour ratio of ~ 1.1 indicates that th e AGN is likely no t domi- 
nating its mid-to-far-IR energetics (see lCoppin et alj20IOl) . 

Given the optical redshifts, the [OI]63 nm line is expected to 
lie between 145.95-146.20 /im. We find a weak ~ 3(t detection 
of [OI]63 centred at 146.51 ± 0.07 /im (corresponding to a red- 
shift of z = 1.319), offset by ~ 500km/s from the Keck optical 
redshift (z = 1.315), which is similar in magnitude to the dif- 
ference between the two zLESS optical redshifts. This shift in the 
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[OI]63 emission from the optical-based redshift is likely a real 
effect since offsets of this magnitude are often seen in BLAGN. 
The best-fitting FWHM of the line emission is 178 ± 113kms"\ 
which is narrower than the instrumental resolution at this wave- 
length (~ 256 km s^^), which is unphysical. We thus proceed to 
fit the spectrum holding the line FWHM fixed to the instrumental 
resolution to yield a less biased line flux (although we note that this 
will artificially decrease the quoted error budget on the fitted inte- 
grated line flux in Table [TJ. In any case, the best-fitting observed 
line width (the line is unresolved) implies that there is no sign of 
an outflow with velocities significantly above the typical Keplerian 
motions of an SMG (~ a few 100 kms~^). This despite the high 
SFR of LESS66 of ~ 370 Mq yr"^ (see Table H w hich w ould 
yield ~ 1 supernova (SN) per year assuming a Salpetea im- 
tial mass function (IMF), resulting in a significant input of mechan- 
ical energy (I SN~ 10** J) into the ISM. 

3.1.5 LESS88 

It has a single secure robust radio counterpart, which is also de- 
tected at 24 fim (formally a tentative counterpart). The ALMA 
data show a single continuum source coincident with this ID. The 
zLESS spectroscopic redshift is z = 1.269 (Q=2). It is not de- 
tected in the X-ray, and has a. 3 a upper limit to its full-band 
flux of So.5-8kcV < 4.1 X 10~^^ erg cm"'^ s~^, derived f rom 
the 250 ks X-ray imaging of the ECDFS jLehmer et alll2005l) . Its 
S8Mm/S4.5^m colour ratio is ~ 0.8, indicating that an AGN is un- 
like ly to be contributing significantly to the energetics of the system 
(see lCoDpin et aLfcoid) . 

Given the optical redshift, the [01] 63 fim line is expected to 
lie at 143.36 fim, and there appears to be some weak emission near 
this position in the spectrum and we claim a tentative detection 
of [OI]63 fim in LESS88. The emission is weak, which precludes 
a successful Gaussian fit. We thus fix the centre and height of a 
Gaussian by eye and proceed to fit for the FWHM (resulting in a 
FWHM of 307 ± 277kms~^ which has been corrected for an in- 
strumental resolution of 261 kms~^ at this wavelength). We note 
that holding the two Gaussian parameters fixed will artificially de- 
crease the quoted error budget on the fitted integrated line flux in 
Table[T] This results in a reasonable line shape and we overplot it on 
Fig- El The line emission is offset approximately by — 300kms^^ 
from the optical redshift, noting that offsets of ~ 200-300 km s~^ 
from the Ha and CO derived redshifts are ofte n seen in infrared- 
bright high-redsh ift star-forming galaxies (e.g. ISturm et"n]|201(]| : 
ISteidel et alj20lol) . 

3.1.6 LESS106 

It has a single secure robust IRAC 8.5 fim counterpart which is also 
detected in the radio (formally only a tentative ID on its own) and 
at 24 fim (not a formal robust or tentative ID). This is the ID we 
targeted with PACS. This counterpart also shows some coincident 
ALMA continuum emission, although we note that the emission is 
rather weak and that the sensitivity of the beam is only ~ 38 per 
cent at this location. There are actually a few weak ALMA contin- 
uum sources present in the beam, and so several sources could be 
responsible for producing the submm flux density seen as a single 
SMG in the LABOCA map. The zLESS spectroscopic redshift of 
our target is z = 1.617 (Q=2). It is detected in the X-ray with a 
full-band flux of Sp s-skcv = (3.4 ± 0.5) x 10"^^ ergcm"^ s"^ 
jLehmer etai]|2005h . which is consistent with t ypical X-ray de- 
tected SMGs (~ 10"^'' - 10"^* ergcm"^ s~^ [Alexander et al.l 



l2005bh . Although there is an AGN in this SMG, its S8;zm/S4.5^m 
colour ratio is ~ 0.8, indicating that an AGN is unlikely to 
be contributing s ignificantly to the energetics of the system (see 
ICoppin et alj20I O). 

Given the optical redshift of our //er.sc/ie/-targeted counter- 
part, z = 1.617, the [OI]63 /im line is expected to lie at 165.34 /xm, 
and a ~ 3 (T emission line is seen at this wavelength (and consis- 
tent with the 1 a error on the optical spectroscopic redshift). The 
line has a FWHM of 413 ± 254kms~^ (corrected for an instru- 
mental resolution of 227kms~^ at the observed wavelength). We 
note that this FW HM is similar to the m ean CO FWHM for SMGs 
(510 ± 80kms~ MBothwell et alj|2012l) . 



3.1.7 LESS114 

It has a single secure robust radio and 24 fim counterpart, which 
was targeted with PACS. The zLESS F0RS2 spectroscopic red- 
shift is z — 1.606 (Q=l). This counterpart also shows some 
ALMA continuum emission. However, there appears to be a nearly- 
equivalently bright ALMA continuum source present in the field, 
and so it is likely that two sources could be responsible for produc- 
ing the submm flux density seen as a single SMG in the LABOCA 
map. It is detected in the X-ray w ith a flux of So.s-g kcv = 
(2.0 ± 0.5) X 10"^^ergcm~2s"^ jLehmer et ai] l2005h . which 
is consistent with t ypical X-ray detected SMGs (~ 10~^® — 
10"^* ergcm"^ s' MAlexander et alj2005bl) . Although there is an 
AGN in this SMG, its S8nm/S4.5^m colour ratio is ~ 0.7, indicat- 
ing that an AGN is unlikely to be contributing significantly to the 
energetics of the system (see lCoppin et alJuOlOl) . 

Given the optical redshift, the [01] 63 fim line is expected to lie 
at 164.65 /im, but there is no apparent emission near this position 
in the spectrum. We therefore derive a 3 cr upper limit to the line 
flux. 



3.2 Far-infrared SED fitting 

We have constrained the far-IR SEDs of our SMG sample by fitting 
their observed submm photometry with a simple modified black- 
body spectrum of the form cx /[exp{h u / k Td) — 1] multi- 
plied by an emissivity function oc u'^, where the dust temperature, 
Td, determines the location of the SED peak and /3 is the power law 
index of our simplified dust emissivity power law. Observations of 
local and distant galaxies constrain /3 to between I and 2 in the 
Rayl eigh- Jeans side of the emission spectrum so we fix /3 = 1. 5 
(e.g.' Hildebriffl3l 19831 : lAgladze et al.lll994lDunne & EdeslbOOlh . 
By fixing (5, we are assuming that the dust properties are similar in 
all the SMGs, but note that we do not account for possible system- 
atics in our quoted errors. For reference, changing the value of /? 
by ±0.5 changes the best-fitting dust temperatures by about ±5- 
10 K and the integrated IR luminosities by ±5-15 per cent. As in 
Section [3TI we ha ve used the l east-s quares minimization MPFIT 
IDL package from lMarkwardtl j2009l) to perform the fitting. 

The submm photometry used in each SED fit includes the de - 
boosted LABOCA 870 /im flux densities from lWeilS et"ai] ( l2009h . 
SPIRE 250, 350, and 500 ^im photometry (see Table |2ll, and also 
our new PACS continuum flux densities or upper limits (from 1 14— 
165 i-im, depending on the redshifted [OI]63 nm location; see Ta- 
ble [TJ. We have extracted the SPIRE 250, 350, and 500 ^im flux 
densities and instrumental errors for our SMGs from the pub- 
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Table 2. We have listed the submm flux densities and our SED fit constraints for the six SMGs in our Herschel-PACS sample. The 250-500 fj,m flux densities 
and instrumental noises are given and hav e been extracted fro m the HerMES SPIRE maps at the SMG counterpart location. We note that we have added 
a confusion noise of 4.8, 5.5, and 6.1 mjy jNguven et alJl20l"ol) in quadrature to the instrumental noise estimat es at 250, 350, an d 500 /^m, respectively, for 
the SED fitting, as per the standard practice. The 870 ^m flux density is the deboosted flux density value from lWeiB et alj j2009l) . The FIR continuum flux, 
LpiR. Tjj and SFR are de rived following the SED fitting p rocedure o utlined in Section \X2\ holding /3 = 1.5 and using the FIR and Lpjp^ definitions from 
[Sanders & Mirabel! il99d) . The SFR is calculated following |KennicutJ 1799^ . 



SMG nickname 


S250 


S350 


S500 




FIR continuum flux 


LpiR 


Td 


SFR 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(xlO"^'' Wm-2) 


(xlOi^Lo) 


(K) 


(Moyr-i) 


LESS21 


26.1 ±3.5 


22.7 ±2.9 


22.4 ±3.8 


7.6 ± 1.3 


1.3 ±0.4 


0.72 ± 0.12 


22 ±2 


100 ± 30 


LESS34 


23.9 ±2.2 


31.0 ±2.1 


26.8 ±2.6 


6.3 ± 1.3 


1.4 ±0.4 


0.59 ±0.09 


21 ±2 


80 ±20 


LESS66 


19.7 ±2.2 


16.1 ± 2.2 


2.9 ±2.5 


5.3 ± 1.5 


4.1 ± 1.2 


1.4 ±0.3 


39 ±5 


370 ± 100 


LESS88 


15.8 ±2.1 


24.2 ± 2.1 


19.9 ±2.5 


4.5 ± 1.4 


1.7 ±0.6 


0.75 ± 0.16 


25 ±3 


140 ± 50 


LESS 106 


24.6 ±2.2 


24.3 ±2.1 


18.1 ± 2.5 


4.0 ±1.4 


2.3 ±0.8 


1.5 ± 0.4 


30 ±3 


350 ± 120 


LESS 114 


42.9 ±2.2 


38.2 ± 2.2 


34.8 ±2.5 


3.9 ±1.4 


5.7 ± 1.0 


3.3 ±0.5 


36 ±3 


850 ± 150 



licly available HerMES DRl mapfl ( iHerMES Collaboration et al.l 
|201^ at the radio or 24 fim counterpart positions. Our SPIRE flux 
density measurements should be less biased on average than other 
methods that could be employed, such as peak matching or measur- 
ing the flux density at the submm position, since there is less com- 
bined uncertainty in the true source position when using the more 
precise r adio or mid-IR positions to measure the SPIRE flux densi- 
ties (e.g. lCoDpin et al. 2008) . We have added confusion noise esti- 
mates of 4.8, 5.5, and 6.1 mJy from lNguven et al ] (l2010h in quadra- 
ture to the instrumental errors at 250, 350 and 500 /^m, respectively, 
and these errors a re used in the SED fit ting routine, as per the stan- 
dard practice (e.g. iMagnelli et alj2012h . 

We then calculate Lpia and the FIR continuum flux using 
the following common d efinitions for ease of comp aring with the 
literature compilation of iGracia-Carpio et alj i201lh who use the 
same definitions. We integrate the best-fitting rest frame SED be- 
tween 40-500 t o derive Lfir, as defined bv lSanders & Mirabell 
( [1996) . Following Sanders & Mirabell ( Il996h . we also calculate the 
FIR continuum flux by integrating the rest frame SED between 
42.5-122.5 fim and dividing this by 47rDL, where Dl is the lu- 
minosity distance. The errors on FIR and Lfir are derived from 
the formal 1 a errors in the fitted parameters computed from the 
covariance matrix returned by MPFIT. Because our single tem- 
perature modified blackbody model does not reproduce the Wien 
side of the far-IR SED well for some of the SMGs (see Fig. [2} we 
are likely to be underestimating the true value of FIR and Lfir 
in these cases . To ch eck this we use the range of SEDs from the 
lDale&Helo3 ( l2002h template library normalised such that they are 
consistent with our measured photometry in order to estimate ap- 
proximate correction factors to our values in Table|2]for each SMG 
(which are of course dependent on the SED templates used). We 
find that our FIR continuum flux and Lfir values for both LESS21 
and LESS34 in Table|2]are potentially underestimated by factors of 
~ 2.8 and 1.7, respectively, and that our FIR continuum flux and 
Lfir values for LESS 1 14 are potentially underestimated by factors 
of ~ 1.6 and 1.5, respectively. Applying these correction factors to 
these SMGs would move these sources down and to the right in 
Fig. [3] by the same factors. We find negligible correction factors 
for our other sources, whose SEDs are better sampled at rest frame 
wavelengths <120 fim (see Fig.|2ll. 

We calculate the SFR for each SMG following iKennicutj 



^ The public SPIRE maps used in this paper were obtained through the 
Herschel Database in Marseille, HeDaM (hedam.oamp.fr/HerMES). 



( Il998l) : SFR(M0yr-^) = 1.7 x 10"'° Lir(L0). This 
relation assumes that the IR luminosity is predominantly 
powered by star formation (i.e. a negligible contribution 
from an AGN, which is a g o od as sum ption in g eneral 
for SMGs; e.g. [Al exander et al.' '2005 a|[3; [Pope et al.l l2008l : 



iMenendez-Delmestre et al. 2009; Coppin et al.ll201(ll ; see also Sec 
tionl3.1 and c omes from a starburst less than 100 Myr old with a 
ISalpeterl ( [l95i) IMF. This relation also assumes that Lir is calcu- 
lated between 8-1000 fim. Because our single temperature mod- 
ified blackbody model does not reproduce the Wien side of the 
far-IR SED nee ded in order to infer the full rest- frame Lir from 
8-1000 u m (e.g. [slain. Barnard & Chapmanll2003h . we use the ap- 
proach of iMagnelli et al. j2012l) for calculating Lir : we have in- 
ferred Lir using the definition o f iHelou et al.l [l988l) , Lfir[42.5- 
122.5 Atm] , and applied a bolometric correction factor of 1.91 from 
[Dale et al.1 (2001). As Magnelh et al. (2012) note, this introduces 
uncertainties in our Lir estimates of the order of ±30 per cent due 
to the variation of intrinsic SED shapes, which is not included in 
our quoted SFR uncertainties. 

Table |2| lists our SED-fit derived FIR, Lfir, Td and SFR 
for each SMG, and Fig. |2| gives our best-fitting SEDs. The SED 
fitting has confirmed that our sample of SMGs are far-IR lumi- 
nous galaxies with Lfir — 0.6-3 xlO'^Lo (and correspond- 
ingly Lir ~ 0.5-5 xlO'^Lo) and Td ~ 20-40 K and harbour 
intense starbursts with SFRs of ~ 100-1000 M© yr"', consistent 
with (while lying towards the lower luminosit y end of the distribu- 
tion) previous SMG SED measurements (e.g. iKovacs et al. [ I2OO6I ; 
ICoppin et alj|2008l ; iMagne fli et al."2012^, and reminiscent of the 
local (U)LIRG population tSanders & Mirabel 1996) . 



3.3 Probing the [OI]63 fim strength in high-redshift SMGs 

It has been an important objective of recent studies to measure the 
strength of the major dense ISM cooling emission lines relative 
to the total FIR emission. Here we investigate this for the impor- 
tant [OI]63 fim '^P2 — Pi fine structure line, which along with 
[01] 146 fim and [CII]158 fim, is a major coolant in the transition 
layers between the molecular and atomic gas in so called PDRs. 
The [OI]63 fim emission line has been relatively unexplored at 
high-redshift, but a well-known 'deficit' or decrease in the strength 
of PDR cooling lines relative to the far-IR continuum emission as 
a function of Lfir has been observed in local infrared-lum inous 
galaxies (e.g. iLuhman et al1l2003l ; [Oracia-Carpio et al]|201 1[) . Our 
new Herschel-PACS data can test whether this deficit extends from 
that seen in local (U)LIRGs to high-redshift SMGs. 



© 0000 RAS, MNRAS 000, 000-000 



Herschel observations of [01] 63 fxm in z 1 SMGs 9 



100 



Q 

X 
3 



a 
Q 

X 



Rest Frame Wavelength [Liml 

100 1000 



LESS21 




100 r 



100 1000 
Observed Wavelength [|j,m] 

Rest Frame Wavelength [um] 

100 1000 



LESS66 




100 1000 
Observed Wavelength [|j,m] 

Rest Frame Wavelength [um] 
100 



1000 



f i00rLESS106 




100 1000 
Observed Wavelength [\im] 



a 

(U 

Q 

X 
3 



Rest Frame Wavelength rum] 

100 1000 




100 1000 
Observed Wavelength [|a,m] 

Rest Frame Wavelength [um] 

100 1000 




100 1000 
Observed Wavelength [|j,m] 

Rest Frame Wavelength [um] 
100 



1000 



100 r 




100 1000 
Observed Wavelength [\im\ 



Figure 2. Best-fitting modified single-temperature blackbody SEDs for our sample of SMGs fit to 250-870 ^m photometry (see Table[2) as well as our new 
PACS continuum flux densities or upper limits (see Table[T). Here [3 has been fixed to 1.5, and the best-fitting and derived values (the FIR continuum flux 
and Lpij^) are listed in Table|2] See Section |X2l for further details. 



The SMGs detected in [OI]63 fj,m emission in our sample have 
an [OI]63 /im-to-FIR ratio of ~0.5-1.5 per cent. Fig. [5] shows the 
[OIJ63 um emission line to FIR continuum ratio plotted against 
Lfir for the SMGs alongside all other low and high-redshift 
galaxy populations with measured [OI]63 j-im emission available 
from the literature. The plot reveals that our SMGs do not show 
an [OI]63 fim deficit, which is different from what is seen in local 
ULIRGs (galaxies with Lfir, ^ 10^^ L0 in Fig. [3j, the purported 
local analogues of SMGs. It thus appears that [OI]63 um is a reli- 
able tracer of star formation in high-redshift SMGs. Our observa- 
tions suggest that local ULIRGs and AGNs are not good analogues 



of typical SMGs, and that SMGs appear to be more like local nor- 
mal/starburst galaxies in terms of their PDR cooling line properties 
relative to their far-IR emission, which we now discuss in detail. 



4 DISCUSSION 

The [CII]158/im fine-structure emission line is the predominant 
coolant of the cold neutral medium, since the C"*" ion is abundant 
and the 158 /im emission line is relatively easy to excite (AE/k ~ 
92 K). In warmer and denser environments (where Go > 10'^ and 
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Figure 3. The observed [OI]63 /jm -to-FIR continuum ratio a s a function of LpiR for SMGs from tliis paper (large black stars) co mpared with avail- 
able Literatu re values (compiled by Gracia-Carpio et al. 2011). The literature compila tion include s all galaxies observed with ISO (Colbert et al. 1993; 
iNegishi et al. 2001 ; M alhotra et a l. 2001 ; Luhman et al. 2003; Dale et al. 2004; Brau her. Dale & Helo u 2008 ) and with PACS ( Gracia-Carpio et al. 20 Ub. 
The hig h-redshift (1.13 < z < 3.04) galaxy comparison sai nple includes 3C 368 iBrauher, Dale & Heloijf2008l) . SMMJ2135-0102 Hvison et al.ll2oiol) . 
SDP.81 jValtchanov et alj201lh . and MIPS J142824.0-H352619 jSturm et al.l2010h . all of which have been corrected for gravitational lensi ng. Several famihar 
well-studied local galaxie s are indicated on the plot for reference, including M82, NGC1068 (a LIRG), and ULIRGs Arp220 and Mrk231 jpischer et al.l2010l ; 
iGracia-Carpio et ^1201 ih . The solid symbols represent the two SMGs with [OI]63 /im detections, and the open symbol represents the SMG which we have 
'tentatively' detected. The open symbols with downward arrows indicate 3 a upper limits. In Section [T2l we found that the FIR continuum flux and Lpjp are 
possibly underestimated by estimated factors of 2.8 and 1.7, respectively, for LESS21 and LESS34, and by factors of 1.6 and 1.5, respectively, for LESS106. 
Applying these coirection factors to these SMGs would move these sources down and to the right in the plot by the same factors. The SMGs do not show an 
[OI]63 /im emission line deficit compared to local ULIRGs, similar to what is seen for the [CII]158 /im emission line for SMGs. 
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n > 3 X 10® cm [OI]63 /^m bec omes the dominant coolant o f 
the interstellar gas (AE/k^ 228 Ki lHollenbach & TielensI 19991) . 
A substantial fraction of the [01] 63 fim emission is believed to orig- 
inate in PDRs at the interfaces between the HII regio ns and molec- 
ular clouds (e.g. iMalhotra et al. 200l| ; IContursi et aL 2002). with 
additional contributions to the line em ission originating from the 
diffuse neutral medium o r shocks (e.g. lSpinoglio & Malkan|[l992l : 
Ivan derWerfetai]! 19931) . [01] 63 fim is thus an ideal tracer of the 
warm, dense, neutral ISM, where star formation occurs. 

The detection of [OI]63 /im in two SMGs in our sample thus 
confirms the presence of a warm dense neutral ISM in these sys- 
tems, and indicates that [OI]63 fim is an important cooling line 
for SMGs. The marginal detection and the non-detections are 
also consistent with this picture. If [OI]63 nm is the dominant 
cooling channel, then the [01] 63 ^m-to-FIR ratio can be inter- 
preted as a proxy of the gas heating efficiency through the dust 
photo-electric effect for a relatively den se gas (n ~ lO^cm"'^; 
[Meijerink, Spaans & Israel 2007; Kaufm an et all 1 19991) . If on the 
other hand [CII]158/^m contributes significantly to the cooling, 
then we can place a lower limit to this heating efficiency. Only such 
dense gas can excite the [OI]63 fim emission line and achieve a rel- 
atively high gas heating efficiency of ^ 1%, since the dense gas 
helps to keep the dust grain charge low and thus the photo-electric 
yield high. 

But what can these observations tell us about the star for- 
mation mode in SMGs, i.e. how their immense luminosities are 
powered? The [OI]63 /im-to-FIR ratios seen in local starbursts are 
~ 0.3%, and even higher ratios are seen i n more 'normal' star 
formation driven galaxies jLord et alJ[T99^ : iMalhotra et a"i]|200ll : 
iGracia-Carpio et ^ 1201 ll) . As mentioned above, local ULIRGs 
tend to show an overall 'deficit' i n this ratio (and al so for other 
important PDR emission lines; e.g. lGracia-Carpio et alj201 ll) rela- 
tive to their less extreme counterparts (see e.g. Fig.[3j. One possible 
cause of this observed 'deficit' is thought to be due to high values 
of the ionization parameter (ie. the number of hydrogen ionizing 
photons, h V > 13.6 eV, per h ydrogen particle) at the surface of the 
clouds (e.g. lAbel et alj2005l ; lAbel et al.ll2009l ; lGracia-Carpio et al.l 
I2OI ih . The effect of increasing the ionization parameter causes a 
larger fraction of the UV photons to be absorbed by dust (which 
are re-emitted in the IR) rather than being available to ionize and 
excite the gas, resulting in a decreased observed [OI]63 /xm-to-FIR 
ratio. Thus, high ionization parameters are the most likely cause of 
the P DR line deficit seen for local (U)LIRGs ( iGracia-Carpio et akl 
l20Ilh . which would indicate that there are intrinsic differences in 
the mechanism or distribution of star formation in these galaxies 
(ie. more concentrated, merger-induced) relative to more 'normal' 
starburst galaxies (ie. more extended, quiescent). 

The fact that we do not see this deficit in high-redshift lu- 
minous star formation dominated galaxies such as SMGs sug- 
gests that the physics dominating the [01] 63 /im-to-FIR ratio 
must be similar to that dominating the star formation activity 
in 'normal' star-forming IR-bright galaxies, which have far-IR 
emission arising from an extended, cooler 'cirrus' component 
(e.g. Efstathiou & Rowan-Robinson 2003). Importantly, the SMGs 
lie about an order of magnitude above the [01] 63 /im-to-FIR ratios 
of local ULIRGs and AGNs, indicating that SMGs are not simple 
scaled up analogues of local ULIRGs, where the star formation oc- 
curs in a compact and highly obscured nuclear region (and poten- 
tially with a significant contribution from an AGN). 

Our result supports what is seen in the [CII]I58 /im emission 
line for SMGs. In a first modest-sized survey of the [CII]I58 /^m 
emission line in a heterogeneous sample of 13 far-IR luminous 



1 < 2: < 2 galaxies, Istacev et al.l ( I2OI0I) find that the luminous 
high-redshift star formation dominated systems in their sample 
(including one gravitationally lensed and one unlensed SMG) do 
not show a [CII] 158 /im deficit; whereas th e AGN-dominated sys- 
tems in their sample do. St acev et alj ( I2OIO ) and De Breuc k et al] 
( [201 l|) use the observed [CII] 158 ^m-to-FIR ratio (which traces 
the far-UV field strength. Go) and scaling argument s to constrain 
the s c ale sizes of the thr ee SMGs to ~ 2-5 kpc l lWolfire et al.l 
'1990). Ivis on et alj j2OI0l) also report a [CII] 158 ^im-to-FIR ratio 
in SMMJ2135-0I02 (a lensed z = 2.3 SMG) which is consis- 
tent with the SM G being powered by s tarburst clumps distributed 
across ^ 2 kpc jSwinbank et"^l2010 t). These studies add to the 
growing body of evidence that high-redshift SMGs are host to ex- 
tensive, galaxy-wide starbursts, in stark contrast to the relatively 
more c onfined starbursts seen in local ULIRGs (typically < 100 s 
pc; e.g. Ilono et al.ll2009l) . These galaxy-wide scale sizes are simi- 
lar to high resolution interferometric mapping of the dense CO gas 
reser voirs of 1 < z < 4 SMGs ( :::::; 4 kpc; Tacconi et al. 20061, 
|2008| ; iBothwelfet alj|20ld ; rEnge] etaHboid ; iHodge et alJboij) 
and of the radio and submm continuum and CO(I-O) emission 



(~ 5 kpc up to 25 kpc; IChaoman et al. 


|2004I Younger et al.ll2008l; 


iBiggs & IvisorJl2008l;lBiggs et al.ll2010 


;Ivisonetal.ll201I). 



Our new [OI]63 pcm measurements have provided a first step 
towards a census of the bright PDR lines in SMGs. We have 
shown that the SMGs in our sample do not show an [OI]63 /^m 
line deficit, similar to what was seen previously in [CII] 158 /im, 
suggestive that this is likely a global property for high-redshift lu- 
minous star-formation-dominated galaxies. These SMGs are thus 
ideal candidates for future observations of CO transitions and of 
the [CII] 158 /im emission line - which together could be used 
to provi de a well-constrain ed PDR solution for the emitting gas 
(see e.g. IStacev et alj [2010 1 to solve for the density, n, and the 
far-UV radiation field strength. Go, in the ISM in SMGs. In- 
terestingly, it is predicted that [01] 63 /im would become more 
important than [CII] 158 /im as the galaxy ISM heats up, which 
is most likely to occur during the most extreme star forma- 
tion events such as in SMGs. At the earliest stages of the star- 
burst event, the C/0 ratio is predicted to evolve strongly - with 
[OI]63 /tm predicted to be strong er than [CII] 158 /im emission 
by at least a fa ctor of ~ 3 (e.g. 'Henkel & Mauersberger 'l993|; 
^leijerink, Spaa ns & Israell i2Q07; Kaufman, Wolfire & Hollenbac^ 
|2007D, since [OI]63 /im traces denser gas (n 3 x 10® cm '') 
than [CII] 158 /im (n ^ 3 x 10^ cm"''). Thus, in principle mea- 
suring the C/O ratio could help to indicate the evolutionary state of 
the star formation episode s of SMGs, using chemical evolutionary 
models such as those of Pi pino & Matt eucci ( 2004). These models 
predict that throughout the starburst episode, the C/O ratio would 
evolve strongly (since C enrichment occurs mainly from low-to- 
intermediate mass stars) and thus this ratio could be used as an 
effective way to 'date' a gala xy and its burst timescale in the early 
throes of star formation (e.g. iMaiolino et al. 2005). Thus, a more 
complete set of diagnostic emission line data would begin to reveal 
a fuller picture of early nucleosynthesis and interstellar processes 
in high-redshift luminous star-forming galaxies, which is critical 
to our understanding of the driving mechanisms involved in these 
modes of star formation. 

5 CONCLUSIONS 

We used Herschel-PACS to target the [OI]63 /im emission line in a 
sample of un/en^ed high-redshift 870 /tm-selected 1.1 < z < 1.6 
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SMGs for the first time. Our sample originally consisted of nine sta- 
tistically robust radio and/or mid-IR counterparts to the SMGs (pro- 
viding precise positions for the SMGs) with high-quality spectro- 
scopic redshifts. We used high-resolution interferometric ALMA 
870 /im continuum data to verify the reliability of our original tar- 
get selection. The ALMA data have verified unambiguously that 
six of the SMG counterparts are bona fide submm-emitting sources; 
three of the SMG counterparts that were targeted here do not appear 
to be a source of submm emission and we have removed these from 
our sample. 

We detect the [OI]63 /im emission line in two of the SMGs, 
tentatively detect [OI]63 emission in one SMG, and measure 
sensitive 3 o upper limits to the line flux in the remaining three 
cases, secure in the knowledge that we have targeted the source 
of the submm emission and the correct redshift. The [01] 63 ^m 
detections confirm the presence of a warm dense neutral ISM in 
these systems and indicate that [OI]63 /im is an important PDR 
cooling line for SMGs. We find that the [OI]63 /im-to-FIR ratio in 
SMGs (ranging from ~0.5-1.5per cent) does not follow the well- 
documented PDR cooling line 'deficit' that has been observed for 
local ULIRGs and some local AGNs which are in a similar far-IR 
luminosity class to SMGs. This result suggests that local ULIRGs 
and AGNs are not good analogues of typical SMGs, and that SMGs 
appear to be more like 'normal' starburst galaxies in terms of their 
star formation properties. This offset in the [01] 63 /xm-to-FIR ra- 
tio compared to local (U)LIRGs has been noted previously for a 
small sample of SMGs in another classical PDR emission line, 
[CII] 158 ^m, and together these results suggest that this could be a 
global property for all the bright PDR cooling lines for SMGs. We 
have taken advantage of a unique window of opportunity provided 
by Herschel to verify that [01] 63 /im is thus a reliable tracer of star 
formation in high-redshift star formation dominated galaxies such 
as SMGs. These Herschel observations are a necessary stepping 
stone to future ALMA studies of [CII] 158 /^m and CO emission in 
order to investigate the ISM physics of SMGs in detail. 
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